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ABSTRACT 

Context. HD 195592 is an 0-type super-giant star, known as a well-established runaway. Recently, a Fermi 7-ray source 
(2FGL J2030. 7-1-4417) with a position compatible with that of HD 195592 has been reported. 

Aims. Our goal is to explore the scenario where HD 195592 is the counterpart of the Fermi 7-ray source modeling the 
non-thermal emission produced in the bowshock of the runaway star. 

Methods. We calculate the spectral energy distribution of the radiation produced in the bowshock of HD 195592 and 
we compare it with Fermi observations of the 2FGL J2030. 7-1-4417. 

Results. We present relativistic particle losses and the resulting radiation of the bowshock of HD 195592 and show that 
the latter is compatible with the detected 7-ray emission. 

Conclusions. We conclude that the Fermi source 2FGL J2030. 7+4417 might be produced, under some energetic assump- 
tions, by inverse Compton up-scattering of photons from the heated dust in the bowshock of the runaway star. HD 
195592 might therefore be the very first object detected belonging to the category of 7-ray emitting runaway massive 
stars, whose existence has been recently predicted. 
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1. Introduction 

'■ The star HD 195592 (DB+43 3630) is a massive runaway 
' visible from the north ern hemisphere, lo c ated at a dis- 
! tance - 1.1 kpc (see ISchilbach fc Roseil (|2008[ )). There 

■ is strong evidence supporting the hypothesis that HD 
195592 is a binary system, w ith a period of about 5 days 

■ (|De Becker. Linder fc Rauwl [2010,). HD 195592 produces 

■ a clearly detected bowshock, as it moves supersonically 
] through the interstellar medium (I SM) as reported by 

■ iNoriega-Crespo. Van Buren fc Dganil (|l997). 

Relativistic particles can be accelerated at strong 
, shocks produced by the stellar wind of a massive run- 

■ away interacting with the ISM (i.e. bo wshocks). These 
' particles can yield non-thermal radiation (jBenaglia et al.l . 
! I20inl : ldel Valle fc Romerol . l2012ll . Recently, a Fermi source, 

2FGL J2030.7H-4417, that might be associated with HD 
195592, has been reported in t he Fermi Large Area 
Telescope Second Source Catalog ([Nolan et al.l . l2012h . see 
Fig. [H In this Lett e r, we apply the model developed by 
Idel Valle fc Romerol ([20T1) for the non-thermal emission 
that takes place in the bowshocks of runaway stars to HD 
195592. We confront then our theoretical spectral energy 
distribution with the measured flux of the Fermi source in 
order to explore the possibility of a physical association of 
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the b owshock associated with HD 195592 (see iPeri et all 
[20TI with 2FGL J2030.7-K4417. 

In Section[2J we make a census of the relevant informa- 
tion related to HD 195592 with emphasis on previous obser- 
vational results in different wavebands. Next, in Sect. [31 we 
discuss the computation of the non-thermal emission and 
present the best-fit to the 2FGL J2030.7-)-4417 source. A 
brief discussion and our conclusions are given in Sect. ID 

2. The stellar system HD 195592 

HD 195592 is an 09.5Ia-type runaway star that is 
thought to be orig inated in the open cluster NGC6913 
([Schilbach fc RoseR [2^008). A detailed spectroscopic study 
revealed that HD 195592 is a binary system with a period 
of a few days, with a lower mass e arly B-type companion 
([De Becker. Linder fc Rauw[ . |2010D . Low amplitude radial 
velocity variations were detected in every strong absorp- 
tion line in the blue spectrum of HD 195592. The variations 
exhibit two time-scales: ~ 5.063 and ~ 20 days. The 5.063- 
day variation is thought to be the period of the binary 
syste m associated with HD 195592. Fo r the second time 
scale, iDe Becker. Linder fc Rauwl ([2OIOI) give two possible 
explanations. It may be the signature of an additional star 
or it may be the signature of intrinsic variability related to 
the stellar rotation. 

Previous radio observations with the Very Large Array 
(VLA) at 4.85, 8.45 and 14.95 GHz revealed thermal emis- 
sion with flux densities of a fraction of mJy, and a spectral 
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Fig. 1. WISE (Wide-field Infrared Survey Explorer) RGB 
image at bands W2 (4.6 /.tm) , W3 (12 ^m) and W4 
(22 /.tm) of the shocked gas around the runaway star HD 
195592; the radiative heating of the swept-up dust pro- 
duces the IR emission, and traces the bowshock. The 95 
% location error circle of the gamma-ray source 2FGL 
J2030.7-f4417 is shown (the 99% countour is outside the 
figure) . The probability countours from gamma-ray sources 
are model— dependent and they must be taken as indicative 
only. 



index equal to .98 (IScuderi et aL. I1998D. Adopting a dis- 
tance of 1.1 kpc (jSchilbach fc Roseil ()2008D ). the integrated 
radio luminosity in the centimeter-domain is estimated to 
be of the order of 4 x 10^® ergs"^. It should be pointed out 
that the radio source reported by IScuderi et al.l ()1998[ ) has 
an angular size of the order of a few arcseconds, and must 
therefore be associated with the stellar system itself and 
not with the more extended bowshock. If the bowshock 
produces radio emission, its flux level should be lower than 
the thermal contribution from the stella r wind s of the 09.5 
and B components, since IScuderi et al.l (|1998D did not re- 



port any extended non-thermal emission. 

In the soft X-ray domain, HD 195592 was never tar- 
get of a dedicated observation. However, the survey for 
point sources in the Cygnus r egion (including t he po sition 
of HD 195592) performed bv iDe Becker et al.l (|2007D with 
the ISGRI instrument on-board the INTEGRAL satellite in 
hard X-rays allowed to derive upper limits for undetected 
point sources. Even though the background level in hard 
X-rays is not uniform, depending notably on the vicinity of 
bright X-ray sources, one could consider upper limits de- 
rived in regions of similar background level to be fairly 
apphcable to the position of HD 19559^ ^- According to 
the flu x upper limit values determined by iDe Becker et al.l 
(|2007D and assuming once again a distance of 1.1 kpc, 
we estimate that our target should not be more luminous 
than about 7 x lO'^^ ergs~^, 5 x lO'^'^ ergs~^ and 7 x lO'^'^ 



^ The upper limits published bv lDe Becker et al.l ('2007') were 
determined in a region located closer to the position of the bright 
X-ray source CygX-3 where the background level should be 
slightly higher than at the position of HD 195592. These val- 
ues should therefore be considered as conservative. This does 
not affect their relevance in the context of this discussion. 



ergs~^, respectively, in the 20-60 keV, 60-100 keV and 100- 
1000 keV energy bands. 

Finally, it might be worth commenting briefly on the 
potential role of binarity in the production of non-thermal 
radiation. It is indeed well established that at least some 
coUiding-wind binaries are able to accelerate particles up 
to relativistic ener gies and co nsequently to p r oduce non- 
thermal radiat i on ( Benaglia fc Romerol 120031 : iDe Beckeil 
[2007; B enaglial |2010D . However, the short orbital period in 
HD 195592 suggests that the stellar separation in the sys- 
tem would not allow relativistic electrons to reach Lorentz 
factors high enough as to produce a significant 7-ray emis- 
sion as detected by Fermi. The particle acceleration process 
would indeed most probably be strongly inhibited by the 
strong ultraviolet/visible radiation fields from both stars 
through inverse Compton (IC) scattering, preventing the 
colliding- winds to emit significantly in the Fermi bandpass. 
In this paper, we will therefore explore the scenario where 
the 7-rays come from the bowshock produced by the stellar 
wind interacting with the ISM. 



3. Non-thermal emission calculation 

For the calculation of the non-therm al radiation we fol- 
low th e model recently developed by Idel Valle fc Romerol 
(j2012t ) . The values for the relevant parameters are given in 
TableHl 

The collision of the supersonic stellar wi nd with the in - 
terstellar medium produces two shocks (e.g. IWilkinl . |2000[ ). 
Relativistic particles are accelerated at the reverse shock, 
that propagates in the opposite direction of the motion of 
the star, inside the stellar wind. This shock is adiabatic 
and strong. The particle acceleration mechanism is diffu- 
sive first order shock acceleration (e.g. iBell Il978l ). The in- 
teractions of the locally injected relativistic particles with 
matter, radiation, and magnetic fields in the shocked wind 
produce non-thermal radiation by a variety of processes 
(Idel Vall e fc Ro merol . l201^ . 

The acceleration region is assumed to be a small region 
near the bowshock apex, of scale length ^ A, where A 
^ M~'^Rq. Here, M is the Mach number of the shocked 
wind and Rq is the so-called standoff radius (e.g. Wilki3, 
" 9961). In the cas e of HD 195592, we adopt Rq - 1.73 pc 



( Peri et al.l . l2012D . 

In order to roughly estimate the magnetic field in the 
fiow, we assume that the magnetic energy density is in 
sub-equipartition with respect to the kinetic energy of 
the winqj. Therefore, we adopt the constraint x < 1- This 
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(1) 



where A is the area of a sphere of radius Rq, and is 
the available power in the system (best fits are provided by 
5 X 10-2). 



X 



The kinetic power of the stellar wind is: 



1 



(2) 



For HD 195592, according to the best available data 
(Table [1]): ^ 10^^ erg s~^. Our estimate relies on the 



Otherwise the gas would be mechanically incompressible. 
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Table 1. Parameters for HD 195592 



Parameter 



value 



Ro 


Standoff radius 


1.73 pc 




Wind mass loss rate 


3.3 xlO"'' Mq yr"^ 


a 


Particle injection index 


2 


T/(a) 


Wind velocity 


2.9x10** cm s"^ 


X 


Subequipartition factor 


5x10"^ 


B 


Magnetic field 


~ 2 X 10"*^ G 




Star temperature 


2.8x10'' K 




Star luminosity 


3.1x10^5 Lq 


Tm 


Dust temperature 


~ 40 K 



Valu es fromlMuiires et all ll2012f) . V alues from 
iMartins. Schaerer fc Hillierl (|2005l ). 



primary star parameters, even though we are dealing with 
a binary system. However, the contribution to the total ki- 
netic power coming from the B-type component is expected 
to amount up to only a small fraction of that from the O 
super-giant, and would therefore not affect seriously our 
order of magnitude estimate. 

The power available to accelerate particles in the reverse 
shock is L = /Lt 2 x 10^'* erg s~^, where / is the ratio 
of the volume of a sphere of radius Rq and the volume of 
the acceleration region. Some fraction of this power goes 
into relativistic particles Lj-d = QrciL. The energetics of the 
gamma-ray source requires Lroi ~ 4 x lO'^"^ erg s~^, so q^ci ~ 
0.2, which seems to be a reasonable value if compared, for 
instance, with supernovae, which are expected to convert 
abo ut 20% of their kinetic power in to relativistic particles 
(e.g. iGinzburg fc Svrovatskiilll964D . 

In the calculations of the spectral energy distribution we 
take into account both hadronic and leptonic content in the 
relativistic power, L^d = Lp + L^- We consider Lp = L^, 
which means equal efficiency in the acceleration of both 
types of particles. 

A proton-dominated scenario seems to be unlikely in 
the case of HD 195592. Since the fraction of the relativistic 
proton energy that goes to neutral pions in each interaction 
is of ~ 17% (e.g. Aharonian & Atoyan, 2000), in order to 
obtain the observed 7-ray luminosity a very high effciency 
in convenrting kinetic energy into relativistic particles is 
necessary. This would required very extended, or perhaps 
even multiple, acceleration sites in the bowsock. For sim- 
plicity, we stick here to the simplest hypothesis: equipati- 
tion beteween electrons and protons. This assumption, from 
the energetic point of view, is also the most conservative. 

The electrons lose energy mainly by IC scattering, 
synchrotron radiation, and relativistic Bremsstrahlung. 
Protons lose energy through proton-proton inelastic colli- 
sions with the shocked wind material. The relativistic parti- 
cles can escape from the acceleration region convected away 
by the stellar wind. These non-radiative losses impose the 
upper limit to the energy of protons. In Fig. [2] we show the 
cooling rates for both electrons and protons in the acceler- 
ation region. The IC scattering of IR photons completely 
dominates the radiative losses. Little power is radiated as 
synchrotron radiation, rendering the non-thermal radio and 
X-ray counterparts quite weak in comparison to the IC 7- 
ray source. 
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Fig. 2. Acceleration and cooling time scales for electrons 
and protons for HD 195592. Up panel is for electrons and 
down panel is for protons. 



As can be seen from Fig. [5] the electrons reach ener- 
gies ^ 0.6 TeV while the protons can get ~ 15 TeV. In 
both cases the Hillas criterion is satisfied, i.e. i?max < 
300(Ai?o/cm)(B/G) eV. 

To determine the steady state particle distributions for 
electrons and pr otons, we solved the transpor t equation 
in steady state (IGinzburg fc Svrovatskil Il964f ). We refer 
to Idel Valle fc Romero! for details of calculation. 

Although most protons are convected away before cooling - 
see Figl5]- downstream instabilities can produce some mix- 
ing with external gas. We consider that a fraction /p (^ 
10%) of the convected protons interact downstream with 
the gas, and hence enhance proton-proton contribution. 
Even if all protons were to cool, the hadronic contribu- 
tion would remain minor, except at energies above 1 TeV. 
At the energies of interest in this paper both the emission 
from the shocked ISM and the absorption are negligible. 

Figure[2] shows the computed spectral energy distribu- 
tion (SED) for the emission produced at the bowshock of 
HD 195592. Fermi data are also shown along with data at 
other wavelengths. T he IR emission from IRAS is produced 
by th e heated dust (jVan Buren. Noriega-Crespo fc Dganl 
[T995h . as the WISE emission. Only the 7-ray flux is non- 
thermal, while at other wavelength the emission is mainly 
thermal. In the calculation we have taken into account only 
the wind of the primary star of the potential binary system, 
so our estimates can be considered as energetically conser- 
vative. 

The SED derived by our calculation lends support 
to the scenario where the 7-ray emission from 2FGL 
J2030. 7-1-4417 could come from the bowshock produced by 
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Fig. 3. Computed SED for HD 195592 bowshock, at 
d - 1.1 kpc and Fermi data of 2FGL J2030. 7+4417. 
The ha rd X-rays upper limits, thermal radio data 
(VL A) (IScuderi et all I1998D and the IR -IRA S- emis- 
sion (jVan Buren. Noriega- Crespo fc Deanil Il995l) are also 
shown. 



HD 195592. Our SED is also compatible with the lack of 
detection of any significant hard X-rays with INTEGRAL 
at the same position. In addition, the predicted synchrotron 
radio flux is too low to have been detected by previous ra- 
dio investigations in the vicinity of HD 195592, where the 
thermal radio emission dominates. 



4. Discussion and conclusions 

Massive stars with strong winds are su spected to be 
7-ray sources smce the early 1980s jC asse fc Paul 



198C 



ay sources smce tne early iij»us llL^asse 
a IVolk fc Formanl . 119821: IChen fc WhitT 



1991 



White fc Chenl. Il 992(1. Despite some statistical evidence 
(iMontmerld . Il979t iRomero. Benagha. fc Torresl I1999D 



conclusive identifications remain elusive to this day. 
This is not surprising taking into account the strong 
non-radiative losses experie nced by rel a tivisti c particles 
in the stellar winds dVolk fc Formanl Il982[ ) and the 
strong absorption expected close t o the massive star 
(|Romero. del Valle fc Orellanal 120101) . The best prospect, 
then, is the detection of high-energy photons wherever 
strong shocks can re-accelerate electrons and ions far 
from the star. This is the case of com bined effects 
of m assive stars in stellar asso ciations ([Torres et aL . 
20041) and colliding wind b inaries (lEichler fc Usovl . 1993 : 



Benaglia fc Romerol l2003t Jbe Becke r, ,200 _; ^. R ecent 
detections of Westerlund 2 f Aharoni an et al.l. 2007) and 
Eta-Carina (jTavani et all . [2009; .Ab do et all . 12010.) . at TeV 
and GeV 7 rays, respectively, seem to support this picture. 

Runaway massive stars offer a unique opportunity to 
detect GeV- TeV emission from single massive stars. The 
stagnation point of the wind of these stars is located 
at sufficient distance as to preclude, under the adequate 
viewing angles, significant 7-ray absorption. The recent 
detection of non-thermal radio emission from the bow- 
shock of BD -H43°3654 by iBenaglia et all (|2010f ) and the 
non-the rmal X-ray emission report ed recently from AE 
Aurigae [L6pez-Santiago et al.l ([20121 ) confirms the capabil- 



ity of some of these stars to accelerate at least electrons up 
to relativistic energies. The presence of rich infrared pho- 
ton fields locally generated by the heated dust swept by the 
shocks guarantees suitable targets for IC interactions, that 
might yield, in some cases, detectable 7-ray fluxes. 

The star HD 195592 presents some characteristics (e.g. 
strong IR field, distant stagnation point due to the rel- 
atively small stellar velocity in a dense medium, absence 
of any other source in the Fermi location error box) that 
makes it a good candidate to be the very first 7-ray emit- 
ting bowshock runaway identified so far. A confirmation of 
the nature of this source would require deep X-ray obser- 
vations to check whether there is a power-law spectrum 
as expected from our modeling. HD 195592 is therefore 
a good candidate for additional observations for instance, 
with ASTRO-H (JAX A mission to be launched in 2014, 
iTakahashi et al.ll2010t ) to investigate non-thermal hard X- 
rays, and with ACIS on the Chandra X-ray Observatory 
because of its low background and high spatial resolution 
necessary to spatially disentangle the soft thermal emission 
from the binary and the expected soft non-thermal X-rays 
from the bowshock. 
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